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ABSTRACT 

The pairing interaction among identical nucleons in a single- 
particle level is treated in the hamiltonian formalism using even 
Grassmann variables. A minimal (irreducible) basis having a re- 
markable symmetry property is set up using composite, commut- 
ing variables with a finite index of nilpotency. Eigenvalues and 
eigenf unctions of given energy, seniority and zero third compo- 
nent of the angular momentum of the pairing hamiltonian are 
then found. The eigenvectors, which cannot be cast solely in 
terms of composite bosons with angular momentum zero and two, 
are expanded in the minimal basis with coefficients analytically 
expressed in terms of a generalized hypergeometric function. 
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A difficulty inherent to nuclear physics relates to the number of nucle- 
ons in a nucleus which is not large enough to exploit the field theoretical 
techniques as successfully as in condensed matter physics, but, on the other 
hand, is large enough to require a high dimensionality for the basis in which 
the physics of the nucleus is treated. 

The standard approach to overcome this obstacle reduces the size of the 
basis by selecting the relevant degrees of freedom, namely those appropriate 
for the description of the low-lying nuclear excitations, and then by introduc- 
ing the associated variables, for example as done in the Arima and Iachello 
model |J. 

Clearly one would like to devise a systematic procedure to carry out 
this reduction starting from an underlying, more fundamental hamiltonian 
expressed in terms of fermionic degrees of freedom. 

Recently we have indeed tried an approach to this problem in the 
path integral framework by performing a non-linear change of variables in 
the Berezin integral defining the partition function of a fermionic system, the 
new integration variables then representing composite bosons of appropriate 
quantum numbers. 

As a simple example, we have applied this method to the problem of 
an even number N of nucleons sitting in a single-particle level with angular 
momentum j (third component m) and interacting through the well-known 
pairing hamiltonian Hp. We have been able to express the pairing action 
and the ground state wave function with zero seniority in terms of the com- 
posites and to recover the familiar formula for the ground state energy of Hp 
[§]. Moreover through this study we learned about several useful properties 
of these composites which have been applied to set up a new perturbative 
expansion in QCD 

However, to solve the problem of the excited states of Hp (those with 
nonzero seniority v ) and to extend the approach to encompass the long range 
quadrupole-quadrupole force as well as the pairing interaction has proven to 
be a hard task to perform. Therefore in this letter we address the « / 
problem without resorting to the change of variables in the Berezin integral. 
Instead we explore whether and how the Fock basis of Hp, set up with 
determinants of single-particle states, can be reduced to a minimal dimension 
and whether the minimal basis can be expressed in terms of composite bosons. 
Indeed one might expect the physics of Hp, which only acts among pairs of 
nucleons coupled to angular momentum J = J z = 0, to be describable by 
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only two composite bosons having the two constituent nucleons coupled to 
J = (the first) and to any allowed non-vanishing angular momentum (the 
second) . 

Specifically, in carrying out the program of reducing the natural Fock 
basis of Hp we shall search for a minimal (irreducible) basis, fully solving the 
pairing problem, using building blocks associated with composite, commuting 
variables having a finite index of nilpotency. Hence these are not to be viewed 
as bona fide bosons, being made up of fermions satisfying the Pauli principle; 
and this is why we shall refer to them as composite bosons. Notably the 
dimensions of the basis will turn out to be fixed by the number n = N/2 of 
pairs present in the problem. 

To illustrate how our approach works, let us start by recalling that in 
the quasi-spin scheme one introduces operators annihilating (creating) pairs 
of particles in orbits having appropriate time-reversal properties, i.e. of the 
type 

[q j 

V m=-j 

(and the hermitian conjugate). In the above (jm,j — m\J0) is the usual 
Clebsch Gordan coefficient and Q = (2j + l)/2. In terms of (0) the pairing 
hamiltonian reads 

H P = -g P BlB Q . (2) 

As mentioned above, one might conjecture that the eigenstates of H P could 
be set up in terms of the operators B and Bj where the specific value of the 
index J(^ 0) in the latter should be irrelevant, since any J ^ corresponds 
to a broken pair. For example, in the spirit of the Arima and Iachello model, 
one could set J = 2. 

However in the framework of the creation and annihilation operators the 
commutator 

CL^ m a m} b} m ib}_ m i = d~ mm > ^1 fl^/flm O^—mO'—m^j j (3) 

which includes, e.g. for J = J' = 0, 

So, Bl\ = n ( 1 - ^ ] , (4) 
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is non-canonical, thus rendering non-trivial the task of finding the eigenstates 
of Hp (n is the fermion number operator). 

An approach that circumvents this difficulty is the hamiltonian framework 
employing Grassmann variables ||. Here one exploits the isomorphism be- 
tween the Fock space F generated by a set of N fermionic creation operators 
a\, ■ • ■ £% and the so-called Q + algebra built with the set of anticommuting 
objects A*, ...X* N (the generators of the algebra). The isomorphism is defined 
by mapping the vectors dj...a]|0 > onto the elements AJ...A*. The image of a 
generic vector |^ >e F under this mapping will be denoted by \1/(A*). Next, 
to a linear operator in normal form, 

6= E E : '---'--' --a i: ..MU,, ...a jk , (5) 

h—ik 3\—3k 

is associated the following function of the Grassmann variables 

0(A*,A)= E O' ■■■■"■^■\] ...K i \ Jx ...\ u . (6) 

h,—ikJi—jk 

which allows one to define the kernel 

K {X*,X) = 0(X*,X) f i + {X*X) . (7) 
Then the action of O on a state \1/ reads 

(OV)(\*) = J [dX'*dX']Ko(X*,X')fi4X'*X')^(X'*) . (8) 

In the above 

At± (A*A) =e ± ^ A * A \ (9) 

In this framework the pairing hamiltonian is then readily expressed in 
terms of the even, nilpotent, commuting Grassmann variables 

Vm = (-iy- m x. m x m (io) 

according to 

3 

Hp = -gpB^Bo with B = fm ■ (H) 

m=l/2 

Now, in seeking a reduction of the basis dimensionality, we shall be guided 
by the two main features of Hp, namely that it 
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1. is expressed solely in terms of the <p, 



2. is invariant for any permutation of the <p. 

These properties suggest also expressing the vectors of the basis in terms 
of the variables (|T0|). For this scope, the action of Hp on the latter should 
be explored. 

This can be accomplished along the lines illustrated in |2[], the result being 

Hpipiy*) = J [d^d^Kpip*,^^^*) = Eip{p*) (12) 

where the integral is over the even elements of the Grassmann algebra as 
defined in [0] and the kernel is 

Kp(p*,r ] )=H P (p*,r ] )e^^ . (13) 

Hence, the above-mentioned difficulty associated with the non-canonical na- 
ture of the commutator (|3]) disappears. 

We then attempt to diagonalize the Hp associated with n = N/2 pairs in 
a basis spanned by states represented as products of n factors p*'s, namely 
of the type 

<Pmi ■ ■ • V*m n ■ (14) 

Since the number of states here is clearly rj , the very large reduction of 
the basis dimension entailed by the choice of the variables (ITO) is apparent: 
indeed in terms of fermionic degrees of freedom, the corresponding basis 
would have a dimension fjn)- 

Moreover, because the variables A anticommute, each vector of the basis 
(|14"D is antisymmetric with respect to the exchange of any pair of fermions, 
and hence the above basis fulfills the Pauli principle. 

We now explore whether, for a given Q and n, the eigenstates of Hp can 
be cast in the form 



m=l 

the index m identifying the set {mi, rri2 • • -?n n } and the j3 m being complex 
coefficients. 
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In the simple situation where n = 1, the states in (ff^) are indeed eigen- 
states of the pairing hamiltonian. In fact the eigenvalue equation, with 
£ = E/g P , is 

(£ + l)p m + E Pp = ° ( 16 ) 

p(^m)=l/2 

and in the basis in (|i~4[) the operator Hp = — ^ + £ is represented by a 
matrix of dimension Q filled by ones but for the principal diagonal, whose 
elements are given by £+1. This matrix is invariant under any permutation 
of the index m which labels its rows and columns. Moreover, of the Q real 
roots of the associated characteristic equation 

{£ + Q) £ n ~ l = 

only two are distinct, namely the lower £o(n = 1) = —Q with multiplicity 
5 = 1 and the upper one £?,{ri = 1) = with multiplicity S = £1 — 1. The 
associated orthogonal eigenvectors are 

*o(n = 1) = 4= E <p* m (17) 

V " m=l/2 

and 

4r 2 („=l)=jV £ (18) 
m=l/2 

(A/" is a normalization factor). Notice that into the eigenvector in (|18|) en- 
ter fi — 1 parameters (the /3 m ), which correspond to the degeneracy of the 
eigenvalue. The above eigenvalues agree with the quasi-spin formalism and 
correspond to the states with seniority v = and v = 2 respectively, v being 
an even non-negative number representing the number of unpaired particles. 

Since the case with n = Q — 1 is identical, but for a shift in energy, to the 
one with n = 1 pairs and, by extension, the case with n pairs is equivalent 
to the one with Q — n pairs, in the following we shall confine ourselves to 
consider only n < ^. 

When n = 2, it is not so trivial to prove that ( |15|) is an eigenfunction of 
the pairing hamiltonian for any Q, although this turns out to be the case. 
For example when Q = 4 the eigenvalue equation is 

(£ + 2)(3 mn + E (P pm + f3 P n)=0 (19) 

p(^n,m)=l/2 



5 



and Tip is represented by the 6x6 matrix whose elements are 8+2 on 
the principal diagonal, on the secondary diagonal and 1 elsewhere. The 
associated characteristic equation 

(£ + 6)(£ + 2) 3 £ 2 = 

has three real distinct roots out of six, namely the lowest So(n = 2) = —6 
with degeneracy 5=1, the intermediate £ 2 {n = 2) = — 2 with degeneracy 
5 = 3 and the highest one £&{n = 2) = with degeneracy 5 = 2. 
The corresponding orthogonal eigenvectors are 

^ ( n = 2) = M (^1/2^3/2 + ^1/2^5/2 + ^1/2^7/2 

+^3/2^5/2 + ^3/2^7/2 + ^5/2^7/2) ( 20 ) 

ty 2 (n = 2) = N 2 [ai (Vl/2^3/2 - ^5/2^7/2) + «2 (Vl/2^5/2 - P3/2V7/2) 

+0 3 (<fy 2 <f*7/ 2 - ^3/2^5/2)] (21) 

# 4 (n = 2) = M [62 ($ya - ^7/2) (^5/2 - ^3/2) 

+^3 (fi/2 ~ <Pt/2) (9*3/2 - f*7/2j} , (22) 

Ao, A2, A/4 being normalization factors and 01,02,03, &2 and 63 free parame- 
ters. 

The case just discussed asks for some comments. First, as for n = 1, 
the matrix representing Tip can be written in Qjl equivalent ways (not all 
distinct), each one corresponding to a different labelling of the states, but 
leading to the same determinant. 

Furthermore in the above example the correspondence between the num- 
ber of parameters that enter into a given wave function and the degeneracy of 
the corresponding eigenvalue is again apparent. In general to an eigenvalue 
£ v (n) of given seniority v correspond || 



(23) 



independent states constructed using the building blocks ([II]) (in our con- 
vention a binomial coefficient with a negative lower index vanishes). The 5 V 
is commonly referred to as the seniority degeneracy. 
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Finally to address the general case we exploit the two lessons learned from 
the previous example. The first one concerns the seniority degeneracy: its 
very existence proves that a further reduction of the dimension of the basis 
can be achieved. The second is that the structure of the symmetric matrix 
of dimension (^j associated with Ttp should be Q 

/ S+n V 1 \ 

(24) 

v V 1 £ + n j 

where the symbol OV 1 indicates that the upper (lower) triangle of the matrix 
is filled with zeros and ones. Indeed the matrix elements oiTip are one when 
the bra and ket differ by the quantum state of one (out of n) pair, otherwise 
they vanish. The diagonal matrix elements simply count the number of pairs. 

An elementary combinatorial analysis then shows that the number of ones 
in each row (column) of the matrix (|24|) is given by n(Q — n). Indeed a non- 
vanishing matrix element has the row specified by n indices whereas, of the 
indices identifying the column, n — 1 should be extracted from those fixing 
the row in all possible ways, which amounts to n possibilities. The missing 
index should then be selected from among the remaining Q — n ones: hence 
the formula n(Q — n) follows. 

To write down explicitly the matrix ( [M[ ) it is convenient to divide the 
set of the Q even Grassmann variables, whose quantum numbers identify the 
levels where the n pairs are to be placed, into two subsets, one with Q — v 
and the other with v elements (to be referred to as I and II, respectively) 
with n < v < Q — n. 

A priori each of these partitions is valid in the sense that it leads to a 
basis with a dimension lower than the one of fllip . On the other hand the 
physics of the pairing hamiltonian is such that for a system of n pairs n + 1 
eigenvalues should be expected, no matter what the value of Q. They of 
course correspond to the breaking of 0, 1, 2 ■ • • n pairs. A partition leading 
to an n + 1 dimensional basis, for any Q and without degeneracy, is the one 
which has v = n. 

Indeed in this instance the rj entries of each row and column of the 
matrix can then be grouped into n + 1 sets, the first one corresponding to 



*As already noticed in the n=2 case, all of the rj ! orderings of the states are equivalent. 
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the n pairs placed in the Q, 
empty (see Fig. la). 



— n levels of I, the remaining n levels of II being 



© 

ii © © 

© © © 



© 

© © 

© © © 

(a) (b) (c) (d) 

Figure 1: The figure shows, in the specific case n=3 and Q=9, the partition 
of the fl levels in the set I (with fl — n=6 levels) and II (with n = 3 levels). 

To this first set correspond do = configurations. To the second set 

are associated configurations with n — 1 pairs in I and one pair in II (see 
Fig. lb), their number being d\ = nr^Zi)- In general the (k + l)-th set 
embodies configurations with n — k pairs in I and k pairs in II, their number 
being 



4 = with < k < n (k integer). (25) 




Clearly ELo 4 =g). 

With this organization of the levels the matrix ( |2^ ) splits into (n + l) 2 
rectangular blocks B^j (with < k,j < n) of dimension dk x dj and, since the 
pairing Hamiltonian only connects states differing by the quantum number 
of one pair, the blocks with \k — j\ > 2 will have vanishing elements. The 
matrix thus becomes block-tridiagonal. 

This suggests the introduction, in lieu of (0), of the following orthonor- 
mal set of n + 1 even commuting variables (the composite "bosons" ) 



\ ^ 



(26) 
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where k varies as in fl2"5|), Sj = J2i=odi, being s_i = and m again identifies 
the set {mi, m 2 • • -m n }. Note that the states (|26|), unlike those in (0), have 
in general an index of nilpotency higher than one. 

The definition fl26|) also reflects our desire to have the composite bosons 
retain as much as possible of the symmetry of Hp. And indeed the 
while not fully symmetric with respect to the interchange of the (p*, turn out 
to be invariant with respect to the interchange of the (p* belonging either 
to set I or to set II. This symmetry property stems from our choice ((26|) 
which enforces the maximum coherence among the components of It is 
remarkable that composite variables corresponding to different combinations 
of the tp* mi if* ■ ■ ■ tp* mn not only hold a lower symmetry than the one displayed 
by (^), but may also lead to the wrong eigenvalues, as we have verified in 
some instances. 

Now in the minimal basis ( |2"E| ) Tip is represented by a (n + 1) x {n + 1) 
matrix whose generic element (M n ) ki =< $£|7^p|$* > results from summing 
the elements of the block B ki of the x matrix, but for the normal- 
ization factor 1 / y/dkdi. The sum is performed by recognizing that all of the 
blocks have the same number of ones in each row. Specifically, in the diago- 
nal block B kk there are (n — k) (Q — 2n + 2k) ones in each row. In the upper 
diagonal block -B^^+i each row instead contains 

c fc = (n - kf (27) 



ones. Since the total number of ones in each row of the matrix (p4|) is n{Vt— n 



the number of ones in the rows of the lower diagonal block B^^-i will be 

b k = k(Q - 2n + k) . (28) 

As a consequence, the non-vanishing elements of the matrix M n turn out to 
be 



(M n ) k;k+1 = < n\n P \n + i >= \j ^-c k , (29) 



fc+i 



and 



(M n ) k+1>k = <<t>l +1 \n P m>=J^b k+1 (30) 



(M n ) kk =< $* k \H P m >=a k = £ + n + n(n-n)-b k -c k . (31) 
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Clearly (M n ) kk+1 = (M n ) k+lk , since the operator Hp is Hermitian and the 
basis (^) is orthonormal. The matrix thus becomes tridiagonal, reading 



/ 



a 

h 




V 



a i 





a-2 



d k-l „ 





d 3 

a k 



c 2 



dk+i K 







\ 



^n— 1 / 



and the associated eigenfunctions should be expanded in terms of the com- 
posite variables (0), namely 



fc=0 A.-0 



(32) 



where the is introduced for convenience. The coefficients of the expan- 
sion (|32"D are then fixed by the eigenvalue equation 



M n u = , (33) 

which, of the equations for the /3's, reduces to n + 1 equations for the 
u>'s. The associated eigenvalues £ obey the secular equation 

D n = det (M n ) = , (34) 

where 8 enters into the diagonal matrix elements a k . 

Now from the general theory of symmetric tridiagonal matrices one knows 
that (|34l) has n+1 distinct and real roots and that these are found by applying 
the recursive relation M 



D n — a n D n _i ; — (c„-i) 2 -D„_2 



dr. 



(35) 



which, after some algebra, yields 



D n = y{y -Q)[y- 2(fi - l)][y - 3(Q - 2)}...[y - n (Q - n + 1)] 



(36) 
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where 

y = n (tt-n)+£ + n. (37) 

We thus see that D n has all of the zeros of D n _\ plus an extra one for 
y = n(Q — n + 1). 

Moreover (|36|) allows us to write down for the general solution of ([34]) the 
expression 

y = p(£l-p+l) withO<p<n. (38) 
Hence the well-known formula for the spectrum of the pairing Hamiltonian 



£ = -(n-p)(p- n -p+l) , (39) 

is recovered, the index p = v/2 being linked to the seniority quantum number 
v. 

For a tridiagonal matrix, a recursive relation among the eigenvectors com- 
ponents, similar to (|3^), can also be established. In the specific case of the 
matrix M n it reads 

d k c k w k+ i = -4-iC/c-iWfc-i - d k a k w k , (40) 

where again < k < n and quantities with negative indices are meant to be 
zero. 

Thus for the lowest eigenvalue y = (p = 0, zero seniority) we have 

(^=0) (v=0) ( v =0) ( A -,\ 

m = w\ = .... = w\ > , (4i) 

namely the collective state 

/Q\ -1/2 n / \ -1/2 (°) 



' fc=0 \ / m=l 



Indeed in (42) all of the components of the wave function, i.e. the monomials 
partition o 



, are coherently summed up. This state obtains for a specific 

n 

the levels defining the matrix M n . However, any other partition 
would lead to the same result, the weights of the components all being equal. 
As a consequence the state ( |42"D has zero degeneracy. 

Concerning the second eigenvalue y — Q (p — 1, seniority 2), according 
to Eq. ([!(]) the components of its eigenstate are 
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w 



(v=2) 



N 2 ( kn - 



(43) 



We have found the following expression for the components of the state 
associated with a generic seniority v 



w 



(v) 



v/2 

^E(-!) j 

3=0 



v/2\ (il-v/2+l 



K 3 f 2 -k, -n + 



V 

2 : 



— n, —n: 1 



(44) 



3_F 2 being a generalized hypergeometric function. In the above the binomial 
(j) is meant to vanish when j > k. 

In particular, for the vector of the maximum seniority w 
sponding to y = n(Q — n + 1) (p = n, seniority 2n), one has 



(v=2n) 



corre- 



(v=2n) 



2/i 



(Q - 2n + k)\(n - k)\ 
{ ' (n-2n)\ 



(45) 



In (|4"3"D, ([Hp and (^) A/2, A4 and A/*2n are normalization constants. 

It is straightforward to show that the general eigenstates (^) reduce to 
the expressions ( |20| , ^1 , ^2|) for Q = 4 and n = 2. 

It is interesting in the present formalism to recover the following impor- 
tant result: for values of Q large with respect to n, the number of dominant 
components in decreases with v, reflecting the weakening of collectivity 
with increasing seniority. Indeed in the limit Q >> n the eigenvalues are 



and moreover 

c k << Q 



y 



2 



Cfc_! ~ kQ and a k 



--k\9. . 



(46) 



(47) 



By inserting the above limits in (f40l) the components corresponding to 
the eigenvalues (|6|) are found to be 



(v) 



w 



(v) 



W 



(v) 

v/2-1 







(48) 
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and 

k 



W k 



l° ] = f^rw^, k = v/2,...,n. (49) 

We thus see that the state with seniority v in the basis of the y/d^^l 
and in the large Q limit has indeed v/2 vanishing components, the remaining 
n—v/2+1 components being expressed through ( fI5| ) via the single component 
In other words, when Q is large the collectivity of a state with seniority 
v decreases as v increases, because its components become fewer and, fur- 
thermore, the surviving components are more capably expressed through a 
single one. 

As a last point in this letter we briefly address the problem of the angular 
momentum of the composite bosons entering in our eigenstates. In this, for 
the sake of illustration we shall confine ourselves to the cases n = 1 and 
n — 2, and the latter only for Q — 4. For this purpose it is natural to 
express the simple building blocks (|TUD, the monomials, as superpositions of 
Q composite bosons Bj of definite even t angular momentum with vanishing 
third component, according to 



(-iy- m ip m = ^Y,{jrn,j - m\J0)Bj . (50) 

Moreover we account for the constraints the nilpotency of the <£>'s variables 
induces on the composite variables Bj. These constraints relate the latter 
among themselves, and hence the Bj are no longer independent. 

It then immediately follows that B£, a composite s-boson, coincides with 
the one pair wave function fllTD . Similarly, with a suitable choice of the 
parameters /3 m , the v — 2 state (|l^) can be made to coincide with £>}, with 
J = 2, 4, • • • 2(Q — 1). This state thus represents a broken pair. 

For two pairs in Q = 4, the v = eigenstate (^) is given by the product 
of two composite s-bosons (two unbroken pairs), namely 

^ (n = 2) = B* 2 . (51) 
Likewise for the v = 2 state (pT]) one has 

ty 2 (n = 2) = (-4ai - 2a 2 + a 3 ) B* B* 2 + 



T Since the ip m are even in m, only even J are allowed. 
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namely a wave function of good angular momentum (J = 2,4,6), as it is 
easily verified by appropriately choosing the parameters. 

Up to this point we have been able to express the eigenstates of Hp only 
in terms of s and d composite bosons. This, however, is no longer the case 
for the v — 4 eigenstate. Indeed, on the one hand it is remarkable that 
Bq cancels out in the differences among the <£>*'s which enter into (|22|), an 
occurrence to be expected since the state v = 4 describes two broken pairs. 
On the other hand, from its explicit expression in term of Bj, namely 



(56 3 - 46 2 ) B* 2 Bl - -j=. (b 3 - 3b 2 ) B* 2 Bl + — y= (lib, + 2b 2 ) B* 4 B; 

(53) 



it appears that (^) cannot be given only in terms of s and d composite 

bosons. Note that the wave functions discussed here are yet to be normalized. 

It is worth pointing out that, since our basis is clearly incomplete insofar 
as the angular momentum is concerned, in general we are bound to obtain 
wave functions that contain superpositions of differing angular momenta. 

In conclusion the aim of this letter has been to illustrate for the specific 
example of the pairing hamiltonian, viewed as a prototype for the interaction 
in complex many-body systems, the suitability of the Grassmann algebra for 
dealing with composite fields. This approach permits a major reduction in 
the number of degrees of freedom used to describe many-body systems, in 
particular to treat the discrete states of nuclei. 

The relevant outcomes of the present analysis arise from the building 
blocks used in setting up our minimal basis. These, in fact, 

• keep track of the nature of their constituents through a finite index of 
nilpotency; 

• while not fully symmetric in the exchange of the <£>'s, nevertheless fulfill 
a remarkable symmetry property. 

Finally, it is seen that more than two composite bosons are required to 
treat the dynamics of Hp. 



* 4 (n = 2) = ~ (5b 3 - Ab 2 ) Bf - A (ll&a + b 2 ) Bf + ^b 2 Bf+ 
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